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Non-equilibrium Relaxation Study of Ferromagnetic Transition 
in Double-Exchange Systems 
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Ferromagnetic transition in double-exchange systems is studied by non- equilibrium relaxation 
technique combined with Monte Carlo calculations. Critical temperature and critical exponents 
are estimated from relaxation of the magnetic moment. The results are consistent with the previ- 
ous Monte Carlo results in thermal equilibrium. The exponents estimated by these independent 
techniques suggest that the universality class of this transition is the same as that of short-range 
interaction models but is different from the mean-field one. 
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Since the proposal by ZenerJ3) the double-exchange 
(DE) mechanism has been studied to explain ferromag- 
netic transition in doped manganites which show the 
colossal magnetoresistance. The ferromagnetism itself 
is well explained by the DE mechanism; at low temper- 
atures, itinerant electrons which interact with localized 
spins through the strong Hund's-rule coupling favor par- 
allel configuration of localized spin in order to gain the 
kinetic energy. The colossal magnetoresistance is under- 
stood by this mechanism also, at least, qualitatively. 

However, in spite of the long history of theoretical 
studies, finite-temperature properties of the DE systems 
have not been fully understood thus far. At finite tem- 
peratures, especially near the ferromagnetic transition, 
fluctuations are quite large through the strong interplay 
between itinerant electrons and localized spins. This 
makes it difficult to describe the thermodynamics quan- 
titatively. Approximational methods, for instance, the 
mean-field approximation, have been known to be insuf- 
ficient because they neglect the fluctuation effects^ 

Recently, the thermodynamics of the DE systems has 
been studied intensively by using the Monte Carlo 
method which fully includes large fluctuations 
A new algorithm is developed by using the moment- 
expansion technique for the density of states to calcu- 
late the MC weight efficiently and the computational 
speed is accelerated considerably by using parallel com- 
puters. This enables us to calculate large-size clusters 
and take account of the finite-size errors systematically 
by the finite-size scaling. The critical temperature T c 
and the critical exponents are estimated by applying the 
finite-size scaling analysis on the MC data. 

The previous MC studies are for thermal equilibrium. 
Some initial MC steps (typically 1,000-10,000 steps) are 
discarded for thermal equilibration, and measurement 
of physical properties is performed by taking the ther- 
mal averages for the grand canonical ensemble using MC 
samples after the thermalization. 

There is another powerful technique to investigate the 



critical properties by using the MC method,4je™the non- 



equilibrium relaxation (NER) technique BLtBlJ' In this 
technique, dependence on MC steps of the order param- 
eter is studied in the non-equilibrium relaxation from an 
initial state which is chosen to be a symmetry-broken 
one. The order parameter decays exponentially above 
the critical temperature T > T c . On the other hand, 
it approaches a constant for T < T c . As a critical re- 
laxation, a power decay is observed at T = T c whose 
exponent relates to the critical exponents. 

One of the advantages of the NER technique is that 
one can study larger-size clusters than by the equilib- 
rium technique since the relaxation is measured in the 
primary part of the MC samplings which are discarded 
as a thermalization process in the equilibrium method. 
The finite-size effect can be exponentially small and ne- 
glected in the relaxation if one calculates a sufficiently 
large-size cluster. The NER technique gives information 
on the critical temperature and the critical exponents, 
independently of the equilibrium MC method. 

In this paper, we investigate the ferromagnetic tran- 
sition in DE systems by the NER technique combined 
with the MC calculations. To our knowledge, this is the 
first application of the technique to critical phenomena 
in itinerant electron systems. The results are compared 
with those by the previous equilibrium technique. 

We study the DE model where itinerant electrons cou- 
ple to localized spins with Ising symmetry on square lat- 
tices. This is a minimum model which shows the fer- 
romagnetic transition at finite temperatures by the DE 
mechanism. We consider the limit of strong Hund's-r 
coupling for simplicity. The Hamiltonian is written b 
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where Cj (cj) annihilates (creates) an electron at site i, 
and 5 describes the localized Ising spin which takes 5 = 
±1. The summation is taken for nearest-neighbor sites. 
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We apply the moment-expansion MC method.© We 
take the same computational conditions as in ref . || . We 
study a square lattice of 24 x 24 sites which seems to be 
large enough to observe the critical relaxation. We will 
comment on the finite-size effect later. The relaxation 
is measured for typically 100-1000 different Markov se- 
quences to estimate statistical errors. The energy unit is 
the half-bandwidth of noninteracting electrons; W = At. 

We calculate the relaxation of the magnetic moment 
per site from the perfectly polarized configuration of lo- 
calized spins. Figure 1 shows the results. For T/W ii 
0.059, the moment decays exponentially. This indicates 
that the system is in the paramagnetic state, i.e., T > T c . 
On the contrary, the data at T/W — 0.056 seem to ap- 
proach a constant, which suggests the broken symme- 
try below T c . Power-decay behavior is observed between 
these temperatures. We estimate T c /W = 0.0575±0.001. 
This result is consistent with the estimate by the equi- 
librium technique, T c /W = 0.058 ± 0.001 B> 
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Fig. 1. Relaxation of the magnetic moment m in Monte Carlo 
steps t from the ferromagnetic state with perfect polarization. 
The data are for T/W = 0.056, 0.057, 0.0575, 0.058, 0.059 and 
0.06 from top to bottom, respectively. The gray lines are the 
least-squares fit by power functions for the asymptotic behavior 
of the data at T/W = 0.057, 0.0575 and 0.058. The dotted line is 
the asymptotic power decay in the mean-field universality class, 



for comparison. 



The asymptotic exponent of the power decay gives in- 
formation on the critical exponents. The magnetic mo- 



The estimate of A is examined to identify the univer- 
sality class of this transition. The mean-field universality 
class gives the exponents as (3 — v = z = 1/2. The dy- 
namical exponent z is obtained from the consideration 
that the MC update from n to n + 1-th step at T = T c 
is described by m n+ i — tanh(m„) in mean-field models. 
Thus, the exponent A is 2 in the mean-field universality 
class. On the other hand, in the universality class of the 
Ising model with short-range interactions in two dimen- 
sions, r-the, critical exponents are estimated as /3 = 1/8, 
v = 1E21E2P and z = 2.16& In this case, the exponent A 
becomes 0.058. Our estimate A = 0.14 ± 0.1 is consis- 
tent with the universality class of the short-range Ising 
model, and is distinct from the mean-field one. 

The present results are compared with the previous 
ones obtained by the equilibrium MC calculations. U 
The previous estimates of the critical exponents were 
/3 = 0.09 ± 0.08 and v = 0.9 ± 0.3, which are consistent 
with those for the short-range Ising model but different 
from the mean-field ones. Therefore, the present result 
by the NER technique is consistent with those by the 
equilibrium technique. It is suggested that the univer- 
sality class of the ferromagnetic transition in the DE sys- 
tem is the same as that of short-range interaction mod- 
els with the same spin symmetry, but different from the 
mean-field one, from these independent techniques. 

We comment on the finite-size effect in the present 
NER study briefly. One of the finite-size effects is an 
exponential decay due to an energy gap inherent in 
finite-size systems. This should decrease the tempera- 
ture where the relaxation shows the power-decay behav- 
ior. Thus, the present estimate of T c is considered to give 
a lower limit. The agreement between the present esti- 
mate and the previous one by the equilibrium technique 
indicates that the finite-size effect is negligibly small in 
the present precision of the data. 

To summarize, we have investigated the ferromag- 
netic transition in the double-exchange systems by the 
non-equilibrium-relaxation Monte-Carlo method. To our 
knowledge, this is the first example of the application 
of this method to the phase transition in itinerant elec- 
tron systems. From the relaxation of the magnetic mo- 
ment, we have estimated the critical temperature and the 
critical exponents. The results are consistent with the 
previous ones which are obtained by the finite-size scal- 
ing analysis on the Monte-Carlo data in thermal equi- 
librium. The exponents estimated by these indepen- 
dent techniques indicate consistently that the universal- 
ity class of this transition appears to be the same as that 
of short-range interaction models but is different from 
the mean-field one. 
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ment m depends on the MC step i as m ~ i in the limit AOYAMA+ project Qhttp: //www.phys.aoyama.ac.jrV 



of t — » oo, where A = [3/vz. Here, (3 and v are the expo 
nents for the order parameter and the correlation length, 
respectively, and z is the dynamical exponent. We fit the 
asymptotic tails of the data at T/W = 0.057, 0.0575 and 
0.058 by power functions as shown by the gray lines in 
Fig. 1 and obtain the estimate A = 0.14 ± 0.1. 
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